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Abstract 

The Brueckner G-matrix for a slab of nuclear matter is analyzed in the singlet l S and triplet 
3 S + 3 D channels. The complete Hilbert space is split into two domains, the model subspace 
5*o, in which the two-particle propagator is calculated explicitly, and the complementary one, 
5", in which the local potential approximation is used. This kind of local approximation was 
previously found to be quite accurate for the 1 S pairing problem. A set of model spaces S (E ) 
with different values of the cut-off energy E is considered, E being the upper limit for the 
single-particle energies of the states belonging to So. The independence of the G-matrix of E 
is assumed as a criterion of validity of the local potential approximation. Such independence 
is obtained within few percent for £"0=10 4- 20 MeV for both the channels under consideration. 
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1. Introduction. 

Most reliable predictions of nuclear properties come from phenomenological approaches 
including the macroscopic-microscopic method |], 0, and, at a more fundamental level, the 
Finite Fermi Systems (FFS) theory PL EL El, the HF method with effective forces |7], H] and the 
new versions of the energy functional method 0. These methods, being comparatively simple, 
permit to carry out systematic calculations for a lot of atomic nuclei. On the other hand, 
being phenomenological, these approaches have the necessity of introducing a set of adjustable 
parameters. This point appears especially delicate when one deals with new types of nuclear 
systems, e.g. nuclei in the drip-line vicinity. Indeed, any extrapolation of phenomenological 
parameters found for stable nuclei is hardly reliable in this case. Therefore the old problem of 
the ab initio calculation of these parameters is of importance not only from an heuristic point 
of view but also from the practical one. 

Application of methods of microscopic theory of nuclear matter to finite nuclei (see, e.g., the 



monographs |nj and [11]) are mainly limited by use of so-called local density approximation 
(LDA). The LDA works reasonably well inside a nucleus but fails at the surface where there is 
a domain of density values for which nuclear matter is unstable. However, it is impossible to 
develop a self-consistent nuclear theory without consistent description of the nuclear surface. 
Indeed, just in the surface region the nuclear mean field sharply changes from zero outside the 
nucleus to a constant value inside it. In terms of the FFS theory, it is associated with a sharp 
variation of the scalar-isoscalar Landau-Migdal amplitude / from a strong attraction outside a 
nucleus to an almost zero inner value. A consistent description of this variation is of primary 
importance for the nuclear theory. We intend to develop an approach to this problem based on 
the Brueckner theory for nonuniform systems beyond the LDA. 

Recently []12[ , dealing with ^-pairing problem for semi-infinite nuclear matter, we developed 
a method of solving the Bethe-Goldstone (BG) type equation for the separable form of the Paris 
A/W-potential without use of any form of local approximation. The same method was applied 



to the case of the slab geometry in Refs. ]T3], H| . Though the separable form of the NN- 
interaction simplifies calculations significantly, they remain to be rather cumbersome and much 
CPU time consuming. To circumvent this difficulty, we devised a new version of the local 
approximation resulting from the investigation of the effective pairing interaction Vg ff [Q. To 
distinguish it from the standard LDA, we named it as the local potential approximation (LPA). 
The effective interaction is associated with splitting of the complete Hilbert space S into two 
domains. The first one is the model subspace So, in which the gap equation is written down 
in terms of the effective interaction Vg ff . The second one is the complementary subspace 5", in 
which the equation for V^ s is obtained in terms of the free AW-interaction V. In this space, 
the pairing effects are not significant, therefore the equation for V^ s has the form of the BG 
equation. Dealing with the pairing problem, the model space was taken in a form convenient for 
nuclear application, which includes all the single-particle states with negative energies e\. The 
LPA is related to calculation of the two-particle propagator in the complementary space which 
enters the equation for V^ s . Roughly speaking, the LPA procedure consists in replacing the 
exact BG propagator by a suitable form taken from infinite nuclear matter. For a fixed value 
of the average centre-of-mass (CM) x-coordinate X = (X 12 + X 34 ) of incoming and outcoming 
nucleons (the x-axis is perpendicular to the layer), the propagator is supposed to be equal to 
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the one of nuclear matter in the potential well V(X). Such an approximation turned out to be 



accurate at a level of a few percent even at the surface \12\. It was shown by a comparison of 



the direct solution for V^ ff with the LPA one. 

From the computational point of view, the problem of finding the Landau-Migdal interaction 
amplitude in terms of the Brueckner G-matrix is much more complicated than the pairing 
problem. First, the additional triplet ^-channel (coupled with the 3 D one) should be considered 
for which the calculations turned out to be more complicated than those for the singlet 1 S'- 
channel. Second, instead of fixing the value of the total perpendicular momentum P± = 
as in the pairing problem, the integral over Pj_ appears in this case. Therefore, though the 



direct numerical solution of the problem, in principle, is possible [15], it looks very difficult. 
Therefore it is natural to attempt using the LPA for solving the BG equation for the G- 
matrix that simplifies calculations significantly. In the case of the singlet ^-channel, the 
accuracy of the LPA for the BG equation is just the same as in the pairing problem because the 
corresponding two-particle propagators in the complementary space coincide. As to the triplet 
channel, applicability of the LPA is not obvious at all. The main goal of this paper is to clarify 
the latter point. 

In the case of the BG equation, there is no evident gain from the introduction of the effective 
interaction and it seems more reasonable, after splitting the Hilbert space as S = So + 5", to 
formulate the LPA procedure in a direct way. According to that splitting, the two-particle 
propagator A in the BG equation can be written as a sum of A = A + A'. The model space 
term, A , should be calculated exactly, whereas the second one, A', within the LPA. It is clear 
that the applicability of the LPA depends on the choice of the model space So- Indeed, all the 
quantum and finite range effects originate mainly from the states nearby the Fermi surface whose 
contribution to the BG equation is strengthened by the small values of the energy denominator 
in the two-particle propagator. These contributions produce the long-range components of A 
and should be taken into account exactly. At the same time, the individual contribution of 
a far- lying state is negligible and only a sum of a big number of such states is important. 
They produce the short-range term of A and can be considered within the local approximation. 
Hence, the accuracy of the LPA should depend on the choice of the model space: wider So is 
higher it is. We use this simple physical idea to impose a criterion of the applicability of the 
LPA. We define the model space So{E ) including all the single particle states with the energies 
E\ < Eq which is more general than that for the pairing problem []. It is obvious that a larger 
Eq corresponds to a higher accuracy of the LPA. We consider the LPA to be valid at some 
value of Eq if the G- matrix does not practically change with additional increase of Eq. 

The paper is organized as follows. Sect. 2 contains the BG equation for the slab of nuclear 
matter with separable iVW-forces. In Sect. 3 the splitting of the Hilbert space on the model 
subspace defined by the cut-off energy E and its complementary subspace is discussed. In 
addition, the LPA for the BG equation is introduced. In Sec. 4 the validity of the LPA for 
determining the G- matrix in the singlet ^-channel is analyzed for different values of Eq. In 
Sec. 5 the analysis is extended to the triplet 3 S + 3 Z)-channel. Sec. 6 contains a summary of 
the results. 



1 The latter corresponds to Eq=0 
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2. Bethe-Goldstone equation for the slab system. 

Let us consider the BG equation for the G-matrix of two nucleons at the Fermi surface, i.e. 
with the single-particle energies e\ = fi, where \x is the chemical potential of the system under 
consideration. In a short notation it reads: 

G(E) = V + VA(E)G(E), (1) 

where V is the free iViV-potential, E = 2/i, and A is the two-particle propagator which is 
determined by the integration over the relative energy of the product (Q P Q P ) of two particle 
single-particle Green's functions. The contribution (Q h Q h ) of two holes is neglected. 

To speed up the convergence it is convenient to renormalize eq. (|1|) in terms of the free 
off-shell T-matrix taken at negative energy E = 2/i. The latter obeys the Lippman-Schwinger 
equation: 

T(E) = V + VA ir (E)T(E), (2) 

where A iv (E) is the propagator of two free nucleons with the total energy E. 
The renormalized BG equation has the form: 

G = T + T(A - A h )G. (3) 

We consider a nuclear-matter slab of thickness 2L placed into the one- dimensional Saxon- 
Woods potential V(x) symmetrical with respect to the point x = with depth of V = 50 MeV, 
the diffuseness parameter d = 0.65 fm, and L = 8 fm. 

We use the separable version \W, [T7| of the Paris iViV-potential [Tj| which for the 1 5'o- 



channel has the (3 x 3) form: 

V(k ) kO=EW* 2 )Ji(* a )- (4) 



For the triplet 3 5*i + 3 .D 1 -channel, a similar (4 x 4) expansion (f|) is valid with a formal 
substitution gi(k 2 ) —>■ (ji(k 2 ) where the column contains two components: 



1, ' (5) 




L being the relative orbital moment in the CM system. 

The scheme of solving the BG equation for a slab of nuclear system in the mixed coordinate- 
momentum representation, which has been devised in Ref.[|15|, is adopted also here. Therefore 



we write down in the explicit form only those equations which are necessary for understanding 
the procedure and refer to []J| for details. We consider first the singlet channel S = 0. As it 



was mentioned above, all relations remain valid for the triplet channel as well after replacing 
9i(k 2 ) — > gi(k 2 ). 

The separable form (^) of the iViV-potential in eqs. (|1|) and (|2]) induces similar expansions 
for the G-matrix 

G(kl, k'l P ± ; Xl , x 2 , x 3 , x 4 ; E) = J2 G^X, X'; E, P±) gi (k 2 ± , x) 9j (k'l, x'), (6) 
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and T-matrix 



T(k 2 ± , k'l, P ± , x u x 2 , x 3 , x 4 ; E) = J2 Tij(X-X'; E, P ±) gi {kl, x) 9j {k'l, x'). (7) 

u 

Here the form factors Qi(k\,x) in the mixed representation are determined by the inverse 
Fourier transformation of gi{kj_+k 2 ) with respect to variable k x . Their analytical form can be 
found in |12] for the singlet channel and in for the triplet one. The obvious notation for the 
CM and relative coordinates in the x-direction are used in (Q) and (0). Of course, the Ty coef- 
ficients depend only on the difference t = X—X' of the CM coordinates. In the perpendicular 
direction, the total momentum and the relative momentum are introduced. 

Substitution of (@) and (^) into eq. (||) results in a set of one-dimensional integral equations: 

G ij (X,X';E,P J _) = T ij (X-X';E,Pi_)+ 
+ dX 1 dX 2 T ll (X -X 1 ;E,P ± )5B lm (X 1 ,X 2 ;E 1 P ± )G mj (X 27 X';E,P ± ) 7 (8) 

Ira 

where 

SBi m = Bi m — Bf m (9) 

is the difference between the convolution B\ m of the two-particle propagator A with two form 
factors gi,g m and the analogous convolution Bf m for the free propagator A iv . The explicit form 
of B im is as follows: 



B lm (X,X';E,P ± )=Y / J 



dk ± (1 -nx) (1 -n\>) 

(2tt) 2 E - Pi/Am -e n - e n > - k 2 Jm X 

xg l nn ,(kl,X)g% n {kl,X'). (10) 



Here we have used a short notation of A = (ri,p_i_), A' = (n',p' ± ), p± = P±/2 + kj_, 
p'_i_ = Pj_/2 — kj_, and n\=9(fi—£\), where E\=e n +p 2 ± /2m, and e n are the eigenenergies of 
the one-dimensional Schrodinger equation with the Saxon- Woods potential. The correspond- 
ing eigenfunctions y n (x) (they are chosen to be real) enter the matrix elements of of the form 
factors: 

sL, n ,{k» ± ,X) = [ dx gi (kl,x)y n {X+x/2)y n ,(X-x/2) (11) 



It should be noted that the symbolic sum over nn' in (|1(]) actually includes the summation over 
discrete states and the integration over the continuum spectrum with the standard substitution 
En^fdp/2n. ^ 

In the singlet channel, the BG equation for the G-matrix is very similar to the one for 
the effective pairing interaction |Tj|. Just as in the latter case, it is convenient to extract the 
singular 5-form Born term from the complete G-matrix: 

G = V + 5G. (12) 

The equation for the correlation component SG of the G-matrix can be readily found from 

SG = VAV + VA5G. (13) 
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An analogous extraction of the Born term should be made also for the T-matrix: 

T = V + ST. (14) 
As a result, the renormalized eq. (|3|) yields the following form: 

5G = F + T{A-A {r )5G, (15) 
where the nonhomogeneous term is 

F = 5T + T(A - A ir )V. (16) 

The explicit transformation of eqs. fll5[) , ( [TBI) to the form similar to eq. fl8|) is quite obvi- 
ous. A simplification of the numerical procedure for solving eq.flTo]) in the slab system under 
consideration can be made using the parity conservation which follows from the symmetry of 
the Hamiltonian under the inversion transformation x — > —x. As a result, the eigenfunctions 
y n can be separated on the even functions, y£, and the odd ones, y~. Then the two-particle 
propagator in the above equations is the sum 

A = A + + A- (17) 

of the even and odd components. The first one, A + , originates from the terms of the sum 
fllpf ) containing states (A, A') with the same parity, and the second one, A~ , from those with 
opposite parity. So long as the AA-potential V does conserve the parity, the propagators A + 
and A~ are not mixed in the BG equation. Therefore the correlation part of the G-matrix is 
also a sum of the even and odd components, 

5G = 5G + + 6G~, (18) 

which obey the separated equations: 

5GT = VA n V + VA n 5G T , (19) 

7T = +,-. 

It is obvious that the integral equation ([H]) can be reduced to the form containing positive 
x values only which simplifies the calculation procedure. This equation should be solved for 
both values of ir separately, then the complete G-matrix is found from eqs. flT2|) , flP8]) . 

All the above general equations remain valid for the triplet channel S — 1. The main change 
occurs in the definition of the convolution integral (|T0|) . For the triplet channel it has the form: 



B?=\X,X';E,P ± )=j:J 



dk ± (1 - n x ) (1 - ny) 

X 



(2tt) 2 E - Pi/Am - £ n - e n > - k 2 Jm 
x ( 9 2(kl, X) gT(kl, X>) + g^(kl, X) g ^{k\, X')) . (20) 

It should be noted also that all multipole expansions of eqs. (|6|) and (|7|) type become the 
2x2 matrices. For example, let us write down the components of the G-matrix in the explicit 
form: 



G LL '(k 2 ± , k'l P ± ; Xl ,x 2 , x 3 , x A ; £)=E Gij (A, X'; E, P ± ) g\ L \^ x)gf'\k'l x'), (21) 

ij 



6 



where L, L' equal to or 2. 



3. Choice of the model space and the local potential approximation. 

The main computational problems of solving the BG equation for the slab system are 
connected with the calculation of the propagators ([T0| ) and (|20|) . The reason for introducing 
the model space Sq(E ), with splitting the complete Hilbert space S — Sq + S' and using the 
LPA in the complementary space S', is as follows. The subspace Sq(Eq) contains all the two- 
particle states (A, A') with both one-particle energies e\,e\> being small, e\,Ey < Eq. In the 
complementary subspace, S'(E ), one of these energies or both of them are large, max(e\, ey) > 
E . For the model space, the contribution of each individual state (A, A') f\to the sum of eq . ([10|) 
or eq.(pO|) is strengthened, in comparison with the analogous one for the complementary space, 
due to a small value of the energy denominator. Such contributions induce the long-range 
terms of the BG propagator A [|15| and should be calculated in a direct way. On the contrary, 



in the complementary subspace no individual state (A, A') is important and only wide intervals 
of the integration over k± contribute to A significantly. The corresponding term of the BG 
propagator is sharply peaked and is mainly determined by the local properties of the system 
T5|. Therefore it is natural to use for it some kind of local approximation. For the problem 



under consideration, it seems to be more natural to use the LPA rather than the LDA because 
the BG propagator in the vicinity of the point X is determined directly by the potential well 
V(X) but not by the density, p{X). At the same time, in the surface region there is no simple 
local relation between p(X) and V(X). 

The splitting of the Hilbert space S = S + S f results in the representation of the BG 
propagator as sum: 

A = A + A', (22) 

where A contains the states (A, A') which belong to the model space, A' including the rest. In 
accordance with the above consideration, we calculate the model space propagator A explicitly, 
but use the LPA for the rest term A'. Obviously, the accuracy of the LPA becomes higher with 
the model space So becoming wider. We consider the LPA to be good at some value of E if 
the calculation results for the G-matrix do not practically change with additional increase of 

Eq. 

The LPA procedure, in principle, is the same for both the channels under consideration 
and is very close to that for the pairing problem [T2L the latter corresponding to the choice 



Eq = 0. Namely, for fixed values of the CM coordinates X12, A34, the convolution integral (|T0| ) 
for S = (or fl2"0|), for S — 1 ) is replaced by the corresponding integral for nuclear matter put 
in the constant potential well Vq = V(X), where X = (X i2 + X 34 )/2, which depends on the 
difference of the CM coordinates t = X u — X 3i : 

B^ A (X 12 ,X 34 ;E,P ± ) = B™(V[X],t;E,P ± ). (23) 

In practice, for a fixed value of the chemical potential \i and the cut-off energy Eq and a 
given set of the potential depths V n , we calculate a basic set Bl^([V n ],t; E = 2/i) of nuclear 
matter propagators. In fact, we used a sequence of V n =SV-(n — 1)) with the step 5V=2 MeV. At 

In fact, the difference e\—[i is small. Just such differences enter the denominator of eq. © at E=2n- 
3 The "individual" state means the fixed value oin,n' and small interval of integration over kj_. 
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a fixed coordinate set Xf., the elements of the LPA propagator matrix Bj^ A (Xi, Xk) are found 
as follows. First, we find the potential depth V(X =(Xi + X^)/2). Then, for a fixed value of 
t = \Xi — Xk\, the LPA propagator is found by a linear extrapolation of two neighboring values 
of £j^([V„], t; E), B'^([V n+ i), t; E), under the condition that the inequality V n < V(X ) < V n+ i 
is satisfied. The convolution integral Bf m for the free propagator A fr , by definition, coincides 
with B\£([Vi = 0], t; E). Details can be found in Refs. Q, p 



4. Validity of the LPA for the singlet channel 

Up to now, we dealt with the general form of the BG equation for the slab system con- 
taining the total perpendicular momentum P± as a parameter. As it was discussed above, 
the "dangerous" terms of the propagator (or of the convolution integrals flUf)) which belong 
to the model space and should be considered explicitly, appear due to the small value of the 
corresponding denominators in the sum of ([K]). It is obvious that they become more dangerous 
if the value of P± becomes smaller. Hence the case of P± = is the most crucial for validity 
of the LPA. Therefore we consider just this particular "bad" case for the LPA. Then, we only 
focus on fi = —8 MeV which is a chemical potential typical of stable nuclei. Thus, we put in 
all the above equations P± — 0, E — — 16 MeV (omitted for brevity from now on). 

As discussed above, the BG equation for the correlation part of the G-matrix (|13|) has a fixed 
parity 7r=(+, — ). Therefore we deal with eq. ([l9|) with fixed 7r which is defined only for positive 
x. After finding the convolution integrals ([H]) and those for the free propagator A iv , the kernel 
of eq.(|l^) and the nonhomogeneous term fllBD are derived by direct integration. Then one 



obtains a set of the integral equations for six independent components of 5Gjj(X, X') (similar 
to eq.(||)) which can be solved numerically [fl2[| , |15|| . Finally, we find the total correlation 



<5G-matrix flT8] ) with components 5Gij(X } X') or, from eq. flT2]) , the complete G-matrix with 
components Gij(X, X'). They differ by a trivial ^-function term: 

5G i3 (X, X') = Gij(X, X') - X i3 5(X - X'). (24) 

The extraction of the latter makes the quantity under consideration more convenient for analysis 
and graphical representation. Therefore, as a rule, we deal with the only correlation part of 
the G-matrix, and not with the complete one. One more remark should be also made before 
discussing the results. Following fli~2"| |, we change the original normalization [^, [T7| of the 



expansion (Q) in such a way that the identity gi(0) = 1 holds true. In this case, the absolute 
values of the A^-coefficients give direct information on the strength of the corresponding terms 
ofthe force. Their values (in MeV-fm 3 ) are as follows: A n = -3.659-10 3 , Ai 2 =2.169- 10 3 , A 22 =- 
1.485-10 3 and Ai3=-2.36-10 1 , A 23 =5.76-10 1 , A 33 =1.72-10 1 . The strengths of all the components 
containing only the indices i = 1, 2 are much stronger (by two orders) than those with the index 
i = 3. Therefore the latter are important only for large momenta which come virtually to the 
BG equation or the Lippman-Schwinger one. If one analyzes the matrix elements of the G- 
matrix over the nuclear wave functions, the typical momenta k ~ kp appear for which the 
contribution of the small components is negligible. Therefore we, as a rule, concentrate on the 
"big" components in a qualitative analysis. Of course, in the calculations all the terms Xik are 
taken into account. 

We made a series of calculations of the G-matrix for several values of the cut-off energy 
£'o = 0, 10, 20 MeV to analyze a dependence of the results on this parameter. To present the 
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results we draw the profile functions 5Gij(X, X'—Xq) of the correlation term of the G-matrix 
at several values of X and the zero moment of the G-matrix: 



G 



(X)= J dtG i:j (X,X+t). (25) 



A typical example of the profile function of 5Gn(X,X'=0) is shown in Fig. 1 for the case 
of the model space with the cut-off energy E =0. It has a sharp peak at the point X—X'. In 
such a scale, similar curves for £"0=10 and 20 MeV are distinguishable from that for E Q =0 only 
after magnification. Such magnified profile functions for large components with ij =11, 12,22 
are shown in Fig. 2, for X'=0, and in Fig. 3, for X'=8 fm. It is easily seen that already 
the difference between the curves for Eq=0 and Eq=10 MeV is rather small. As to that for 
£ =i0 MeV and E =20 MeV, it looks negligible. 

To analyze the dependence of the G-matrix on Eq in a more quantitative way, it is worth 
to compare each other the zero moments ( 25|) at different values of Eq. They are shown in 



Fig. 4 for the same large components and, as an example, for one small component, ij=13. 
One sees that a difference, at a level of a few percent, exists between the curves for E =0 
and E =10 MeV and again the additional increase of Eq from 10 MeV to 20 MeV does not 
practically influences the results. This is true not only for big components, but also for small 
ones. 

Finally, we calculated the "Fermi averaged" G-matrix in the ^-channel: 

< G F > 5 =o (X) = J £G ij {X)g i (4(X))g j {4(X)), (26) 

ij 



where we have introduced the local Fermi momentum as fcp(X)=y 2m(//— V(X)) at fi—V(X) > 
and which otherwise takes zero value. Such an average appears if one calculates the Landau- 



Migdal amplitude in terms of the G-matrix [|15 |. To this respect, one remark should be made. 



Though the profile functions Gij(X,X') are strongly peaked at the point X = X', the long 
range "tails", which are hardly seen "by eyes" in Fig. 1, also contribute to the zero moment 
(|25|) .f\ These terms of the G-matrix appear due to states entering the model space and their 



contribution to the integral (25) was analyzed in [T5| for the case of Eq = 0. When one deals 



with the problem of evaluation of the Landau-Migdal amplitude which is supposed to be a short- 



range coordinate function it is natural to cut these tales. In Ref . |n| a recipe was suggested to 
use the Fermi averaged G-matrix by eq.(^) with the zero moments "with cut-off" which are 
defined by the integral of the (|25|) type, but with limits of \t\ < t c , t c =3 fm. Of course, for the 
problem of the validity of the LPA it is not important what kind of the zero moment is used in 



(|26|). However, we use here the same recipe for the Fermi averaged G-matrix as in because 
it is more physical. This quantity is shown in Fig. 5 for the same three values of Eq together 
with the analogous average value of the free off-shell T-matrix: 

< Tp > s=0 (X) = Y t f ij (E=2^g i {t4{X))g j (4(X)), (27) 



lr This contribution depends on ij and is, as a rule, not greater than 10 -j- 20 %. 
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where the zero moments Ty of the T-matrix are defined similar to (^J). In this case, the 
introduction of the cut-off with t c =3 fm does not practically change the integral. Of course, it 
is X-independent. 

Again, the difference between the Fermi-averaged G-matrix for E o =10 MeV and that for 
E =20 MeV is negligible. Their deviation from the one corresponding to E =0 is also very 
small everywhere except in the surface region. It should be noted that the difference between 
the average G-matrix and T-matrix is rather small. A similar property was found previously 
13| , for the effective pairing interaction in the ^-channel. 



Analysis of Figs. 2-5 leads us to the conclusion that for the singlet channel S = the LPA 
works perfectly well for E Q =10 4- 20 MeV. Moreover, within the accuracy of a few percent, it is 



also valid for E =0. The latter agrees with the analysis of |L2| where the LPA was introduced 
for the pairing problem in the ^-channel. 

5. Validity of the LPA for the triplet channel 

In general, the calculation scheme for the triplet 3 S + 3 T>-channel is very similar to that for 
the singlet one, though the calculations become more cumbersome in this case. Indeed, first, 
we get ten independent components Gij(X, X') and ten integral equations (|Sp for them instead 
of the six of the singlet case. Second, the calculation of the convolution integral fl2"0|) in the 
triplet channel is also more difficult than that of fllQf) . Therefore the problem of simplifying 
these calculations is even more important than in the singlet channel. 

Contrarily to the singlet case, now it is difficult to separate the multipole terms into the 
"large" and "small" ones. Again we changed the original normalization [|TjJ of the expansion 
(|),(|D to guarantee the identity #f =0 (0) = 1 (it should be noted that #f =2 (0) = 0). Then the 
strengths of the corresponding terms of the force (in MeV-fm 3 ) are as follows: An= — 1.618 x 
10 3 , Ai 2 = - 1.296 x 10 3 , Ai 3 =8.921 x 10 2 , A 14 =4.271 x 10 1 , A 22 =7.848 x 10 2 , A 23 =1.394 x 
10 3 , A 24 = - 7.860 x 10 2 , A 33 = - 7.450 x 10 2 , A 34 = - 5.723 x 10 2 , A 44 =1.865 x 10 3 . These 
values show that, though the strengths of different components vary significantly, only one of 
them, Ai 4 , is less by two orders of magnitude as compared to the largest ones. Therefore almost 
all the terms are significant. We take several typical components to illustrate the calculation 
results. 

The profile functions and zero moments are shown in Figs. 6-8. One can see that now the 
results with increase of the cut-off energy from E =0 to E =10 MeV change more sizably than in 
the singlet channel, especially in the surface region. At the same time, the subsequent increase 
of Eq up to 20 MeV does not practically influence the G-matrix, the maximum variation being 
of a few percent. Hence, once more one may conclude that the LPA is sufficiently accurate if 
the cut-off energy is E =10 4- 20 MeV. But, contrarily to the singlet case, the accuracy of the 
LPA is rather poor when the model space is limited to E =0, since at the surface the G-matrix 
must tend to the off-shell free T-matrix. But the latter has a virtual pole at small energy. It is 
then clear that an accurate account of the contribution of the single-particle states with small 
positive energies is important for correct description of this pole behavior. Therefore these 
states should be included into the model space 5*0. This does occur if one chooses the cut-off 
energy E > 10 MeV, but it does not occur if one takes E =0. 

Let us now consider the Fermi averaged G-matrix in the triplet channel, which now is a 
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2x2 matrix in the orbital angular momentum space: 

< G F (X) = Y,G!=\X)g^{kl{X))gf\kl{X)l (28) 

where L, L'—O, 2. Just as in the singlet case, the quantity Gjj in fl28|) has the meaning of the 
zero moment "with cut-off". The components of this matrix are shown in Fig. 9 for all three 
values of the cut-off energy Eq. The component < Gp >°s=i is significantly larger than those 

(2) 

containing L=2, especially at the surface region where the form factors g\ vanish. Again all 
the components of the Fermi-averaged G-matrix calculated for E =10 MeV coincide practically 
with those for £"0=20 MeV, though deviations from the £"0=0 case can be noticeable. 

6. Conclusion 

The applicability of the LPA has been analyzed for the Brueckner G-matrix. Previously 
12 1 this kind of the local approximation proved to be quite accurate in the problem of the 
microscopic evaluation of the effective pairing interaction in the ^-channel. The BG equation 
for a slab of nuclear matter has been solved for the singlet 1 S and triplet 3 S + 3 D channels using 
the separable representation |L6], |17[ of the Paris potential. The complete Hilbert space has 
been split into two domains separated by the energy Eq. The model subspace Sq(Eq), in which 
the two-particle BG propagator is calculated explicitly, contains all the two-particle states with 
both the single-particle energies £\,£\< < E . In the complementary subspace, S'(E ), the LPA 
for the BG propagator has been used. A qualitative analysis shows that the accuracy of the 
LPA becomes higher with increase of the energy Eq. It should be also higher for larger values 
of the perpendicular total momentum P±, therefore we limit ourselves to the most "dangerous" 
case of P±=0. 

For either channel under consideration, a set of calculations of the G-matrix has been made 
for different values of the cut-off energy Eq. The LPA has been assumed to be valid starting 
from the value of Eq for which the G-matrix does not practically change any longer. An 
approximate independence of results on the value of Eq, at a level of a few percent, was found 
for £q=10 -T- 20 MeV for both the channels. It should be mentioned that in the singlet channel 



the accuracy of the LPA is sufficiently high even at E =0, in accordance with |[L2| . On the 
contrary, in the triplet channel the LPA is not practically applicable at E =0. 

A similar analysis could be made also for the channels with L > 0. Estimates show that 
in this case conditions for validity of the LPA are even better than those for L=0. We believe 
that the LPA gives a good device for microscopic describing finite nuclear systems. 

In this paper, we limit ourselves to one value of the chemical potential /i= — 8 MeV which 
is typical for stable nuclei. In principle, for smaller values of \x the analysis should be repeated 
However, as some estimates show, even in the drip-line vicinity where /x — > the LPA should 
be rather good for Eq=10 -j- 20 MeV in either channel. At the same time, at Eq=0 it should 
become inapplicable even in the singlet channel. 
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Figure captions 

1. The profile function 5G U (X,X' = 0) in the singlet channel for E = 20 MeV. 

2. The profile functions 5Gij(X,X' = 0) in the singlet channel for E = (dotted lines), 
E = 10 MeV (dashed lines), and E = 20 MeV (solid lines). 

3. The profile functions 5Gij(X,X' = 8) in the singlet channel for E = (dotted lines), 
E = 10 MeV (dashed lines), and E = 20 MeV (solid lines). 

4. The zero moments Gij(X) in the singlet channel for Eq = (dotted lines), Eq = 10 MeV 
(dashed lines), and E = 20 MeV (solid lines). 

5. The Fermi-averaged G-matrix in the singlet channel < G-p >s=o {X) for E = (dotted 
line), E = 10 MeV (dashed line), and E = 20 MeV (solid line) and the Fermi-averaged 
T-matrix (thin solid line). 

6. The same as in Fig. 2 for the triplet channel. 

7. The same as in Fig. 3 for the triplet channel. 

8. The same as in Fig. 4 for the triplet channel. 

9. The same as in Fig. 5 for the triplet channel. 
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